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ABSTRACT 



We analytically and numerically investigate the possibility that a still undis- 
covered body X, moving along an unbound hyperbolic path from outside the 
solar system, may penetrate its inner regions in the next few years posing a 
threat to the Earth. By conservatively using as initial position the lower bounds 
on the present-day distance dx of X dynamically inferred from the gravitational 
perturbations induced by it on the orbital motions of the planets of the solar 
system, both the analyses show that, in order to reach the Earth's orbit in the 
next 2 yr, X should move at a highly unrealistic speed v, whatever its mass Mx 
is. For example, by assuming for it a solar (Mx =M ) or brown dwarf mass 
(Mx = 80mj up ), now at not less than dx = 11 — 6 kau (1 kau=1000 astronomical 
units), v would be of the order of 6 — 10% and 3 — 5% of the speed of light 
c, respectively. By assuming larger present-day distances for X, on the basis of 
the lacking of direct observational evidences of electromagnetic origin for it, its 
speed would be even higher. Instead, the fastest solitary massive objects known 
so far, like hypervelocity stars (HVSs) and supernova remnants (SRs), travel at 
v ~ 0.002 — 0.005c, having acquired so huge velocities in some of the most violent 
astrophysical phenomena like interactions with supermassive galactic black holes 
and supernova explosions. It turns out that the orbit of the Earth would not be 
macroscopically altered by a close (0.2 au) passage of such an ultrafast body X 
in the next 2 yr. On the contrary, our planet would be hurled into the space if 
a Sun-sized body X would encounter it by moving at v/c = 10 -4 . On the other 
hand, this would imply that such a X should be now at just 20 — 30 au, contrary 
to all direct observational and indirect dynamical evidences. 

Subject headings: gravitation; celestial mechanics; planet-star interactions; methods: 
analytical; method: numerical 
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Introduction 



Several free-floating astronomical bodies traveling in the interstellar space in the Milky 
Way have been recently detected. 

In recent years a handful (16) of unbound astrophysical objects lonely wandering 



light, have been discovered ( 


Brown et al. 


2005: 


Edelmann et al. 


2005; 


Hirsch et al. 


2005; 


Brown et al. 


2006a 


b 




2007a 


b; 


Heber et al. 


2008; 


Brown et al. ' 


2009; 


Tillich et all 20091: 


Brown et al. 


2010: 


Irrgang et al. 


2010). Thev are the so-called hvpervelocitv stars (HVSs 



whose existenc e as a consequenc e of the Massive Black Hole (MBH) host ed in t he cen ter 



of the Galaxy flGhez et all 120081 : iGillessen et all 120091 ). was predicted bv iHilld ( 119881 ). 



Gravitational mechanisms of ejection based on three-body mutual interactio ns of binary 
syst ems with the MBH, or possibly a pair of MBHs, have been proposed by iHillsl ([1988 ) 



and lYu fc Tremaind (120031) . The consequent rates of HVSs cr eation would be of the order 
of 10" 3 - 10" 4 yr" 1 flPerets et al 1120071: lYu fc Tremainel [20031). About 10 3 HVSs may exist 
within the Galactic solar circle ( lYu fc Tremainel 120031 ) . Contrary to those neutron stars 
exhibiting high proper motions, which are supernova remnants (SR s), known HVSs are 
mostly B-type main-sequence stars. As an example, HE 0437 — 5439 (lEdelmann et al.ll2005l). 
movi ng at heliocentric speed v = 723 km s~ l = 152.517 au y r" 1 = 0.0024c (lEdelmann et al. 
20051) . is a B star with mass M = 9M (IBrown et al.l l2Q10h . The study by iBrown et al. 



( 120101 ) has yielded the first compelling evidence that these HVSs actually come from the 
center of the Galaxy. All the known HVSs are at about 50 kpc and are unbound with 
respect to the Galaxy. 

Another class of isolated astrophysical objects moving at very fast speeds, not related 
to HVSs, is represented by those neutron sta rs which are t he remnants of asymmetric 
explosions of core-collapse supernovae (S Ne) (IBurrowsll2000l ). Their extreme speeds are 
very likely to be attributed to the kick ( Fryerl 2004 ) received in such a kind of peculiar 
deflagration^]. By measuring the displacements of young pulsars from the apparent 
centers of their associate d SN shells and using the pulsar spin-down periods as age 
estimates, ICaraveol (119931 ) and iFrail et al.l ( 119941 ) inferred that pulsars are typically born 



with transverse velocities of 500 km s -1 , and that velocities v > 2000 km s _1 may be 
possible. At present, the observational record belongs to the radio-quiet neutron star RX 
J0822 — 4300, which moves at a record speed of 1570 km s _1 = 331.191 au yr _1 = 0.0052c 
at a distance of 7000 lyr = 4 x 10 8 au, as measured in 2007 by the Chandra X-ray 
Observatory (jWinkler fc Petrell2007l ). It is thought to have been produced in an asymmetric 
SN explosion. 



Hndeed, if the explosion of a progenitor star expels the ejecta preferentially in one direc- 
tion, the compact core must recoil in the opposite direction because of momentum conser- 
vation. 
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Moving to isolated substellar objects having smaller velocities by about one order of 
magnitude (v ~ 10 _4 c), we have the so-called brown dwarfs. They are astrophysical objects 
in the range mass M ~ 0.04 — 0.09M o = 41 — 94mj up unable to sustain hydrogen fusion in 
their cores; as a consequence, it is very difficult to detect them, since most of the energy of 
gravitational contraction is radiated away within 10 8 yr, leaving only a ve ry low r e sidua l 
lum inosity. After that their e xistence was postulated for the first time by iKumarl (119631 ) 
and lHayashi &: Nakanol (119631). the first undi sputed discovered brown dwarf, and the first 
T dwarf, was Gl 229B ( Nakajima et al.l 1995 ). with a mass M = 20 — 50mj up . After the 
advent of large-area surveys with nea r-infrared (IR) ca pability in the late 1990's, hundreds 
more brown dwarfs were discovered (iKirkpatrick 2005 ). Actually, s maller brown dwarf, 



In particular, in 2005 iLuhman et al 



with M < 20m Jup , exist (ILodieu et al.l 120071 ) 
discovered Cha 110913 — 773444. It is a planetary-mass brown dwarf with M 



(120051 ) 

8m Jup , 

which is well within the mass range observed for bounded extrasolar planets (M < 15m j up 
An even smal ler body, named rh o Oph 4450 with M = 2 — 3mj up , has been recently 
discovered by iMarsh et al.l (120101 ) . 



Concerning the existence of free-floating planets of smaller mass, IStevensonl (119991 ) 
noted that, under certain circumstances, Earth-sized solid bodies wandering in the 
interstellar space after being ejected during the formation of their parent stellar systems 
may sustai n forms of life. Agai n as a consequence of three-body interactions with Jovian 
gas giants, iDebes fc Sigurdssonl (120071 ) have recently shown that during planet formation a 
non-negligible fraction of terrestrial-sized planets with lunar-sized com panions will likely be 
ejecte d into interstellar space with the companion bound to the planet. IDebes fc Sigurdsson 
(120071 ) yield a total number of free-floating binary planets in the Galaxy as large as 7 x 10 s . 
At present, no planets like them have yet been detected. Proposed microlensing surveys o f 
next generation will be sensitive to free-floating terrestrial planets (IBennett fc Rhid 120021 ) ; 
under certain circum stances, they may be able to yield 10100 detections of Earth-mass 
free-floating planet s (IBennett fc Rhid 120021). One to a few detections could be made with 
all-sky IR surveys (IDebes fc Sigurdssonl 120071 ) . 



Are there some solitary traveling astronomical objects, still undetected for some 
reasons, which may hit the Earth over a time scale of a few years? In view of the growing 
attention that such a possibility may really occur orj^l 21 December 2012 is receiving in 
larger portions, also (relatively) educated, of the large public, the present study may also 
have a somewhat pedagogical/educational value contributing, hopefully, to dissipate certain 
fears too often artificially induced simply for the sake of gain. Mere academic disdain 
and/or conceit, derision, and hurling insults should not be retained as adequate practices to 
counter them. Moreover, the analysis presented here can be repeated in future when other 
"doomsday" dates will likely pop out. 



2 See, e.g., 



ht t p : / /en . wikipedia . org/ wiki / Nibiru_collision 



on the WEB. 
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The paper is organized as follows. In Section [2] we present a relatively simplified 
analytical calculation^ which, however, grasp the essential features of the situation 
investigated. A more sophisticated numerical analysis is presented in Section |3j It is based 
on the numerical integration of the equations of motion by randomly varying the initial 
conditions. Section [4] summarizes our findings. 



2. Analytical calculation 



Let us consider a simplified two-body scenario in which a test particle X moves along 
a heliocentric hyperbola hurling itself towards th e Earth. Its conserved (positive) total 
mechanical energy E is ( ILandau fc Lifshitzlll976[ ) 

= 2^ ~ r > °' (1) 

where r and v are the relative X-Sun distance and speed, respectively, \x is the system's 
reduced mass 

• M X M 

f 1 = iT7 , —* » ( 2 ) 



M x + M 



and 



a = GM X M G 



0) 



(3) 

where G is the Newtonian constant of gravitat ion. The semi-major axis a of the hyperbola 
is determined by its total energy according to ( ILandau fc Lifshita 119761 ) 



a 

2E' 



(4) 



The eccentricity e > 1, which, in general, depends on E and on the conserved orbital angular 
momentum L, can be fixed by making the simplifying assumption that the perihelion 
distance of X 

r (p c ") = a (e - 1) (5) 
coincides with, say, the average heliocentric distance of the Earth 

„2 



(r e ) = a e ( 1 + -f ) = 1.000142 au. 



Thus, 



1. 



(6) 



(7) 



3 See also |http: //www.badastronomy.com/bad/misc/planetx/orbitmath.html on the 
WEB for the case of a bound, highly eccentric orbit of X coming close to the Earth. 
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Thus, e depen ds now only on the conse rved energy E. The parametric equations for the 
hyperbola are (ILandau fc Lifshitzlll976l ) 



t 



a (e cosh£ — 1) , 
^(esinhe-0 



(8) 



where the parameter £ takes all values from — oo to +oo; at perihelion £ = 0. Let us, now, 
fix r to a given value. It is the heliocentric distance dx at which the putative X should 
be located at the present epoch. Remaining in the realm of celestial mechanics, dx can 
be thought as dynamically constrained by its perturbations of the orbital motions of the 
known bound major bodies of the solar system. In particular, upper limits on the tidal 
parameter of X 

Kx = -jr- (9) 



have been recently obtained (jlorioll2009l ) by using the secular precessions of the longitudes 
of the perihelia w of the inner planets: for each assumed value of the X's mass there is 
a different lower limit for dx- Although, strictly speaking, they have been obtained by 
assuming X fixed during a planetary orbital revolution, we will use them for the sake of 
concreteness. Of course, if we quite reasonably postulate that X is made of baryonicQ matter 
emitting electromagnetic radiation, other, tighter bounds on its present-day distance may 
be derived from its electromagneti c direct detectabilit y. The recently launched Wide-field 
Infrared Survey Explorer (WISE) (IWright et al.ll2010l ) will survey the en tire sky in the 
mid-IR with far greater sensitivity than any previous all-sky IR surve ys (IPricd 120091) l ike, 
e.g., that performed by the Infrared Astronomical Satellitqj (IRAS) ( Beichmanl 1 198 71 ) . 
Among the scientific goals of WISE there is also the detection of solitary brown dwarf-like 
bodies in the neighborhood of the solar system. WISE should be able to reveal the existence 
of a body with the mass of Jupiter withiro 63 kau = 0.3 pc = 1 lyr, while a lightweight brown 
dwarf with Mx = 2 — 3mj up would be detectable up to 412 — 618 kau = 2 — 3pc = 7— 10 lyr. 
Moreover, WISE could find a Neptune-sized object out to 700 au. Now, by keeping r = d x 
it is possible to extract the contemporary value of the parameter £ corresponding to the 



4 Concerning the putative existence of stars and planets made of a particu lar ki nd of 



non-ba ryo nic dark matter, i. e. th e so-called mirror matter ( )Lee fc Yangl Il956l ). see Foot 
f ll999aM : Foot fc Silaeadzel (l200lh . 



See http://spider.ipac.caltech.edu/staff/tchester/iras/no_tenth_planet_yet.html on the 



WEB about the alleged discovery of a planet in the remote peripheries of the solar sys- 
tem a by IRAS. 

6 See on the WEB http://www.planetary.org/blog/article/00002070/. 
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present-day distance of X 



£ = arccosh 



dx 
a 



(10) 



By substituting eq. f lTUj) into the parametric equation of t of eq. ([8]), one can plot the time 
required to pass from dx to (r e ) as a function of the alleged velocity of X at the present 
epoch: it is sufficient to evaluate E of eq. ([T]) for r = dx- Thus, from the value of the 
velocity required to take a given time interval-typically of a few yr-to reach the Earth's 
orbit starting from dx, it is possible to make reasonable guesses about the plausibility of 
the hypothesis that such a putative body X moving towards our planet actually exists out 
there. 

To be more specific , let us assume that X is an object with the mass of the Sun; 
in this case, lloriol (120091 ) yields dx = 12 kau as dynamically inferred lower bound of its 
present heliocentric distance. Figure [1] shows that such a Sun-sized body X should travel 
at a implausibly high speed (v ~ 0.06 — 0.1c) to reach our orbit in a few years from now. 
Recall that the highest recorded speeds of unbound objects of stellar size are 0.002 — 0.005c. 

M X = M Sun 



3 
2.8 



V 2.6 



2.4 
2.2 
2 



0.06 



0.07 



0.08 



0.09 



0.1 



Fig. 1. — Time t, in yr, required to a body X with Mx =M Q to reach the terrestrial orbit 
from dx = 12 kau as a function of its present day speed fix, in units of c. 



Note that the situation is even worse if we take a larger value for the limit distance dx in 
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view of the fact that, after all, a baryonic star should have been easily detected if it was 
really at just 12 kau from us. Indeed, it turns out that by setting, say, dx = 100 kau the 
required speed would closely approach c. It may be of interest to note that, by traveling at 
v = 0.002 — 0.005c, a Sun-sized body X would take 300 — 800 yr to reach our orbit if it 
was now at 100 kau from us, while 40 — 90 yr would be required if it was at just 12 kau. 
Incidentally, let us remark that the closest black hole so far discovered, whose distance has 
been directly measured from its parallax using astrom etric VLBI observations, is in the 
X-ray binary V404 Cyg, at ab out 2 kpc = 4x 10 8 au ( jMiller- Jones et al.ll2009l ). Another 
close black hole is V4641 Sgr ( lOrosz et al.l 120011 ). at about 7 — 12 kpc. 



Figure |2] depicts the case of a brown dwarf with M = 80mj up and dx = 5.2 kau ( llorio 

20091 ). Also such a scenario looks highly implausible because it should be v/c m 0.03 — 0.05; 



M x = 80m Jup 




0.03 



0.035 



0.04 

fix 



0.045 



0.05 



Fig. 2. — Time t, in yr, required to a body X with Mx = 80mj up to reach the terrestrial 
orbit from dx = 5.7 kau as a function of its present day speed (3x, in units of c. 



moreover, for dx = 20 kau it turns out that v/c ~ 0.08 — 0.2. No brown dwarfs at all 
moving at speeds comparable to those of SNRs and HV Ss are known; on the contrary, their 
speeds are of the order of v ~ 100 km s -1 = 3 x 10 _4 c f lFaherty et al.ll2009l ). Traveling at 
such typical speeds, it would take 1 — 3 kyr to reach the terrestrial orbit for dx = 20 kau, 
and 300 — 900 yr for dx = 5.7 kau. 
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Figure [3] shows that also the case in which the putative colliding X is a rock-ice bo dy 
with the mass of the Earth is unlikely because, by assuming dx = 175 au (Ilorioll2009l ). 
it should travel at v/c = 0.001 — 0.003 to reach the orbit of our planet in the next few 
yr. Given the ejection mechanisms occurring in the planet formation processes which may 



M x = m E 



ar 




0.001 



0.0015 



0.002 
fix 



0.0025 



0.003 



Fig. 3. — Time t, in yr, required to a body X with Mx = m e to reach the terrestrial orbit 
from dx = 175 clU. clS db function of its present day speed flx, in units of c. 

be responsible for such free-floating small planets, their typical velocities should be of 
the order of v ~ 1 — 3 km s _1 = 1 — 0.3 x 10 _5 c for a Jupiter-sized mass ejecting body 

300 yr would be required by traveling at such speeds if 
175 au. 



flGoldreich et al.ll2004h . Thus, 180 
an Earth-sized body X was now at dx 



As we will see in Section [3j the conclusions of such a simplistic analytical two-body 
scenario are also supported by a more sophisticated, numerical analysis. 

It may be interesting to note that some reflections by M. Brown similar to the reasonings 



developed in detail in this Section can be found at |http: / /news .discovery. com/space/mike-brown- planetx-pfi 
on the Internet. The case of a body, of unspecified mass, reaching the Earth's orbit on an 
unbound trajectory in the next 2 yr starting now from 1 kau is touched. Strictly speaking, 
the speed of such an unbound X is computed by as suming that it travels uniformly, so that 
it is v = 2.4 x 10 3 km s _1 = 0.008c. According to lloriol (120091 ). 1 kau is the dynamically 
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inferred lower limit for a body with Mx = mj up lying perpendicularly to the ecliptic; the 
speed required to come here in the next 2 — 1.6 yr turns out to be 0. 8 — 1% of c. If we 
take dx = 1-2 kau for a jovian-sized body lurking now in the ecliptic (jloriol 120091 ). we get 
v/c = 0.01 to reach 1 au in the next 2 yr. Concerning a Jupiter-sized body X, Brown 
at 



http://news.discovery.com/space/mike-brown-planetx-pluto.html puts it at at a few 

2.5 kau we have v/c = 0.02. 



thousand au; in this case, by setting, say, dx 



3. Numerical calculation 

We, first, numerically integrated the equations of motion of an unbound body X in a 
ICRF/J2000.0 heliocentric frame with a coordinate system employing rectangular Cartesian 
coordinates along with the ecliptic and mean equinox of reference epoch J2000. In regard to 
the initial position chosen for X, we took the predicted coordinates of the Earth at to =21 
December 2012 retrieved from the HORIZONS WEB interface by NASA/JPL and added 
randomly generated small corrections to them, i.e., 

Xq ^ = x ffi (to) + 5 X1 

Zq X — z ®(to) + &zi 



where 6 x ,5 y ,5 z were randomly generated from a uniform distribution within ±0.001 au. 
Concerning the initial velocity, we randomly generated it by imposing the conditions 



v (X) > / 2G(M S +M X ) 



^ \l 

(12) 

«o W < c, 



where 



2G(Me + A« (13) 

is the limit parabolic velocity; a hyperbola occurs if v > v p . Starting from such sets of 
randomly generated initial conditions, we numerically propagated the trajectory of X 
backward in time over 2 yr, so that ta n represents the present-day epoch. In such a way, 
by performing several runs, the conclusions of Section [2] turned out to be substantially 
confirmed in the sense that, in order to avoid finding X at the end of the integration, i.e. at 
the present epoch, closer than the dynamically inferred lower limits dx, too high velocities 
would be required. 
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Then, we made a further numerical analysis in which we used the final state vectors of 
X of each of the previous runs backward in time as initial conditions for new runs performed, 
now, forward in time over 2 yr. In other words, now to corresponds to the present epoch, 
while tfi n = 21 December 2012. In such new runs we also added the Earth, Jupiter and 
Saturn by modeling their mutual interactions and their attractions on X. Their initial 
conditions, corresponding to the present epoch, were retrieved from the HORIZONS WEB 
interface. The situation remains unaltered: starting today from positions corresponding 
to the dynamically inferred lower limits dx, all the numerically propagated trajectories of 
X reach heliocentric distances of about 1 au in next 2 yr traveling at unrealistically high 
speeds, as seen in Section |2j It turns out that, also according to such an analysis, larger 
initial distances for X yield even larger speeds for it, just as in Section [2j The inclusion 
of the major planets of the solar system do not cause noticeable alterations to such a 
picture. Conversely, from our numerical analysis it turns out that the hypothetical passage 
of such a fast body X would not distort the orbits of the planets considered, in particular 
of the Earth. This is clearly depicted by Figure |4] which shows the numerically integrated 
terrestrial orbit in the next 2 yr in the case of a Sun-sized X body supposed located today 
at 11.241 kau and moving with v/c = 0.076. Incidentally, in the example showed the 
mutual X-Earth distance at tfi n amounts to Ar = 0.2 au. Instead, a much smaller velocity 
of X would induce macroscopically noticeable changes in the Earth's orbit, as shown in 
Figure El It is obtained for v/c = 2.2 x 10~ 4 , with v/v p = 1.1. Such a scenario would be 
catastrophic since in it the Earth would be finally stripped from its orbit and thrown away, 
as an extension of the time span of the numerical integration to 5 yr shows. Of course, it is 
highly unrealistic since it implies the present existence of an undiscovered Sun-sized body 
X at just 26.3 au. 



4. Summary and conclusions 

We analytically and numerically investigated the possibility-which cyclically gains 
popularity for a variety of psychological and/or sociological reasons in extended portions 
of the large public, even cultivated-that a yet undiscovered astronomical body X, moving 
on an unbound trajectory from outside the solar system, may penetrate its inner regions 
by closely encountering the Earth in the next few years. For the sake of concreteness we 
choose a time span of 2 yr ending at 21 December 2012, familiar to a non-negligible amount 
of people, but the strategy outlined here can naturally be extended to any temporal interval 
and dates in the not unlikely case that in the more or less near future-presumably after 
2012-other analogous "doomsdays" of astronomical origin will be proposed. 

As initial positions, we conservatively choose the lower limits dx for the present-day 
distance of such a putative X from the bounds dynamically inferred from the magnitude 
of the perturbations that it would induce on the orbital motions of the inner planets of 
the solar system. Given that, at present, there are no direct observational evidences of 
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Fig. 4. — Section in the {xy} plane of the numerically integrated Earth's orbit over the next 
2 yr by assuming for X Mx =M , dx = 11.241 kau, v = 0.076c. 

electromagnetic origin for the existence of X, tighter constraints on its distance, i.e. larger 
values for dx, may well have been adopted. The initial velocities were chosen by allowing 
for unbound, hyperbolic trajectories in the field of the Sun. Both analytical and numerical 
calculations, performed for different values of the mass of X by randomly varying its initial 
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Earth's orbit: {xy} section 



1 




-1 -0.5 0.5 1 

x (au) 



Fig. 5. — Section in the {xy} plane of the numerically integrated Earth's orbit over the next 
2 yr by assuming for X Mx =M , dx = 26.3 au, v = 2.2 x 10~ 4 c. 

conditions and including also the Earth, Jupiter and Saturn, show that, in all cases, X 
should move at unrealistically high velocities to reach heliocentric distances of 1 au in the 
next 2 yr. No known astrophysical objects with high speeds, acquired in certain known 
physical processes, move as fast as the putative X should do. In the case of a body with 
the mass of the Sun or of a typical brown dwarf (Mx = 80mj up ) the speed required to 
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come close the Earth in the next 2 yr from presently assumed distances of thousands-ten 
thousands astronomical units would be 6 — 10% and 3 — 5% of c, respectively. Even higher 
speeds are involved if we adopt larger values of the initial distance of X relying upon the 
still missing direct detection of it from electromagnetic radiation. The fastest Sun-sized 
objects known so far travel at speeds as large as 0.2 — 0.5% of c, and are produced in some of 
the most violent astrophysical processes known like interactions with supermassive galactic 
black holes and supernova deflagrations. Moreover, it turns out that the orbit of the Earth 
would not be distorted in a macroscopically noticeable way by the close (0.2 au) passage of 
such a hypothetical ultrafast body. The terrestrial path would be sensibly altered in such 
a way that the Earth would be thrown away if the speeds involved by a passing star-sized 
body were quite smaller, of the order of 0.01% of c. Such a scenario is highly unrealistic 
because, in this case, X should be now at just a few ten astronomical units. 
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